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Introduction 



The existense of a huge gradient of the gravitational 
field in massive astrophysical objects leads to radial sep- 
aration of matter . In other words, in a strong gravita- 
tional field there are strata with the different equations 
of state. Each strata has its equation of state. In par- 
ticular a neutron star has a structure of this kind. 

A relativistic fluid ball with an inhomogeneous static 
stratified matter configuration is considered. The ball 
is filled with a Pascal perfect fluid. The distribution of 
matter in the star is supposed to be spherically symmet- 
ric. The processes connected with rotation and radia- 
tion are excluded. The interior model with a stratified 
structure of mass density is constructed. 

We use the geometric system of units with the veloc- 
ity of light c = I and the Newton gravitational constant 
Gn = 1, i.e. x = 87r. The metric is taken in Bondi's 
form 

ds 2 = G{r) 2 dt 2 + 2L(r)dtdr 



-r 2 (d9 2 + sin 9 2 d<p 2 ), 



(1) 



where G{r) 2 and L(r) are functions from the radial 
variable r; t is the time coordinate; and tp are angle 
variables. Further we will denote d/dr as a prime. 
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2. The Einstein equations 

The gravitational field, described by the metric , can 
be found from the Einstein equations 



R *-l R 9* = -"' T «» 



(1) 



where Rik is the Ricci tensor and R is the scalar cur- 
vature. 

The energy-momentum tensor of a Pascal perfect 
fluid can be written as 



T ik = (p (r) +fJ.(r))- UiU k - p (r) g ik , 



(3) 



where p = p(r) is the pressure, fj, = fi(r) is the density 
of matter, and it, is the 4-velocity. 

The gravitational equations in dimensionless vari- 
ables can be written in the form convenient for numerical 
calculations after elementary transformations: 



e(x) = I ■ / a(x) ■ x 2 dx; 
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(4) 



•G = 0; (5) 



\ • ■'• • 1> + —7- ) ' (M + P) > ( 6 ) 



where x — r/R is the dimensionless radius; differentia- 
tion in x is denoted by a prime; R is the radius of the 
astrophysical object; \ — K 'R 2 \ 



e{x)=G 2 {x)/L 2 (x). 



(7) 



In our case, the following matching conditions at x 
1 can be written 



e{x = 1) = i - 77; G{x = l) = v/T^?; 



(8) 
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G'(x = l)= ; (9) 

p(x = l) = 0, (10) 

where r\ = 2m/ R is the compactness factor, and m is 
the stellar mass. 

3. Stratified structure 

The four functions Fix) , L(x) , p{x) , and (J,(x) are un- 
knowns, and we must specify one of them. For example, 
we take /j,(x). 

We will add the mass density profile to Einstein's 
equations. This function can be constructed as a multi- 
step continuous function. The mass density distribution 
is simulated here by a function of this kind ([1], [2]) 

r \ /„ sin(a • x) 2 \ 

H(x) =/io- ( l-b-x + — -^— L - h (1) 

where /jq is the central density, a is a dimensionless 
parameter (0 < a < oo) which controls the number of 
layers (see Fig.l), & is a dimensionless parameter con- 
nected with the stellar properties. 
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Figure 1. The mass density distribution. There are 
two strata for a = 5 and three strata for a — 10, the 
parameter b = 1 is fixed. 

The functions G{x) , L(x), p(x), K(x) =p(x)//j,(x) 
are found numerically by the Runge-Kutta-Fehlberg 
method. Here we must require the validity of the dom- 
inant energy condition, \p(x)\ < p{x), p > 0, /i > 0. 
Matching to the Schwarzschild exterior solution is car- 
ried out at j] = 0.1 . 

These functions can be plotted for different values of 
the parameter a with fixed values of b. For example, 
we present a plot of p(x) (see Fig. 2). 
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Figure 2. The pressure. The parameter values are a — 
5, o = 10, 6= 1. 

4. Neutron star 

Now we will consider a neutron star model and find the 
compactness r\ for such an object. 

We have the following equation for parameters a , b 
from the condition (8) 

\\{x)x*dx = % (12) 
Jo xR 

Its solution depends on the stellar mass m and ra- 
dius R, and the central density po . This mass density 
is connected with the energy density /io by /io = po • c 2 , 
where c is the speed of light; here and henceforth we 
work in the CGS units. 

We take the characteristic physical parameters of a 
neutron star: 

R = 1.5 ■ 10 6 cm, m = 1.9891 • 10 33 g, (13) 

po = 4.8064 • 10 1 Vcm 3 . (14) 

Then the solution of Eq. (12) gives a = 10, for b = 1. 

The star is filled with neutrons. Accordingly we will 
use the parametric equation of state of degenerate fermi- 
gas: 

A) = (sinh(0-0 K/c 2 ; 
Po = (sinh (0 - 8 sinh (£/2) + 3£) K/3, (15) 

where K = rrtfC 5 I (327r 2 fi, 3 ) , nif is a rest mass of 
fcrmion, h is the Plank constant, and £ is a parame- 
ter, < £ < oo . 
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Let us find the mass density /j,q and the pressure po 
at the stellar center from this equation. We have for 

po = 1-0038 • 10 li g/cm 3 , 

Po = 1.07 • W 33 dyn/cm 2 . (16) 

The compactness factor is found from Eq.(12) for the 
parameter values a = 10 and /zo = 1-0038 • 10 14 • c 2 . As 
a result, rj w 0.16. On the other hand, the characteris- 
tic compactness factor for neutron stars found from the 
observations is rj w 0.1. 

5. Gravitational phase transitions 

Besides, we must note that such a model with strat- 
ified structure has the Petrov algebraic type changing 
from layer to layer. In other words, an extremum point 
exists in each layer, with the first and the second deriva- 
tives of the mass density are equal to zero. In particu- 
lar, such a point is the centre of the ball under study. In 
neighborhoods of such points, the metric function g 00 
will coincide with that in the interior Schwarzschild so- 
lution with a high precision, with its own integration 
constants. It means that the gravitational field belongs 
to the algebraic type (a conformally flat type). In 
the other points of the model, the gravitational field be- 
longs to the algebraic type D. Thus in the interior part 
of our model the algebraic type of the gravitational field 
varies: into the centre (and its neighborhood) is the type 
which passes on to type D. 

From the point of view of phase transition theory, 
we have a second-order phase transition. The phases of 
a "substance" here are the algebraic types of gravita- 
tional fields [3]. On the other hand, we have an cusp 
catastrophe from the point of view of the catastrophe 
theory. Further on, in the the second layer there is also 
a second-order phase transition in the gravitational field 
from type D to type and again to type Z>,and so on 
up to the ball surface. The Schwarzschid exterior solu- 
tion belongs to the type D as we know. 

In other words, as we move from the center to the 
surface, an alternation of the algebraic types of gravi- 
tational fields by the second-order phase transitions is 
observed. 

Thus the properties of the model confirm the theo- 
rem from [4] which says that any spherical gravitational 
field can only belong to two algebraic types: D and 0. 

Another example of phase transitions of the gravi- 
tational fields as transitions from one algebraic type to 
another can be found in [5]. 

6. Summary 

In conclusion, we must note that the model of an as- 
trophysical object (a neutron star ) with a stratified 
structure has been constructed as a result of studing 



a relativistic liquid ball. This model allows for obtain- 
ing a qualitative description of real stars with stratified 
structure. Furthermore, in such stars we can expect the 
existence of the gravitational phase transitions as the 
alternation of Petrov's algebraic types. 
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